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Abstract 

This paper deals with the location of the complex zeros of the Tutte 
polynomial for a class of self-dual graphs. For this class of graphs, as the 
form of the eigenvalues is known, the regions of the complex plane can 
be focused on the sets where there is only one dominant eigenvalue in 
particular containing the positive half plane. Thus, in these regions, the 
analyticity of the pressure can be derived easily. Next, some examples of 
graphs with their Tutte polynomial having a few number of eigenvalues 
are given. The cases of the strip of triangles with a double edge, the 
wheel and the cycle with an edge having a high order of multiplicity are 
presented. In particular, for this last example, we remark that the well 
known conjecture of Chen et al. [2] is false in the finite case. 

1 Introduction 

The program of studying complex zeros of the partition function was pioneered 
by Yang and Lee [7] . The location of the complex zeros of the Tutte polynomial 
is important from a statistical mechanic point of view because it is related to a 
possible phase transition for q-Potts model on this family of graphs using the 
Fortuin Kasteleyn representation (for some correlation duality relations for the 
planar Potts model, see [!)). There exists an impressive literature concerning 
the location of zeros of the Tutte and chromatic polynomials for a large class 
of graphs and, in this context, the Beraha numbers play an important role (see 
for example [5] and references therein) . Sokal [6] proves for general graphs that 
these zeros are dense in the complex plane. A large number of conjectures for 
a wide class of graphs such as triangulation, planar, cubic graphs are given in 
[3]. For self dual graphs, an intriguing conjecture proposed by Chen and al. [2\ 
asserts that, in a half plane, the complex zeros are located on a circle. 

First, our goal was to prove this conjecture for the large class of self dual 
graphs as possible. In fact, in the finite case, we believed that the self duality, 
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implying the symmetry of the Tutte polynomial, was the adequate property in 
order to obtain such a result. More precisely, we thought that it forced in the 
positive half plane the complex zeros of this polynomial to be located on the unit 
circle in the same spirit as the Lee and Yang theorem. We realized that it is not 
true. We give an example and surprisingly we find a sequence of self dual graphs 
for which this unit circle does not belong to the accumulation set of the zeros. 
But the question raised by this conjecture remains at the thermodynamic limit. 
Unfortunately in the general case, for example for self dual strips of the square 
lattice, the eigenvalues coming from the powerful transfer matrix method are 
roots of polynomials with high degree: then, it is difficult to study the location 
of curves of degeneration of the dominant eigenvalue even if we focus on the 
positive half plane. 

In this work, we choose to study a family of self dual graphs. For this class 
of graphs, as the form of the eigenvalues is known, we are able to present regions 
having only one dominant eigenvalue in particular containing the positive half 
plane. Then we derive easily in these regions of the complex plane the analyt- 
icity of the pressure. Next, we study some examples containing a few number 
of eigenvalues. We use deletion and contraction rules of the Tutte polynomial 
to obtain recursive formula rather than the transfer-matrix method. One orig- 
inality of this approach is to unify in the same framework different kind of self 
dual corrections. Besides, we use a suitable variable allowing us to identify the 
set of accumulation of zeros more easily. 

The outline of this paper is the following. In the first section, the conjecture 
of [2] and the well-known results such as Beraha Kahane Weiss theorem and 
Vitali's convergence theorem are recalled. Next, general results, giving regions 
where we have only one dominant eigenvalue, are established and the analyticity 
of the pressure is deduced. Then for the strip of triangles with a double edge, 
the wheel, the cycle with an edge having a high order of multiplicity, the set 
of degeneration of the dominant eigenvalue is obtained. To conclude, some 
perspectives and conjectures raised by ours results are proposed. The technical 
aspects of the proofs are given in the last section. 

2 Preliminaries and framework 

In this paper, we deal with self dual graphs (G = G*) where G* denotes the 
dual graph of G. One useful tool in graph theory is the Tutte or dichromatic 
polynomial. It is a function of two variables x and y; we denote T{G,x,y) the 
Tutte polynomial of the graph G calculated at the point (x,y). A well known 
fact is that for self dual graphs, this function is symmetric 

T(G,x,y) = T(G*,x,y)=T(G,y,x). 

On a statistical mechanic point of view, the link with Ising and g-Potts models 
using the Fortuin Kasteleyn representation leads to look in particular on the 
hyperbola (x — l)(y — 1) = q. In all the following, we assume that q > 1. 
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Usually, we consider the variable v — re l9 by taking 

x = 1 + and y = 1 + v / 9 w - 

In our case, we choose working with the variable 

z = x + y — 2 = + 1/u) then = z + q + 1 . 

It is easy to see that, on this curve and from the symmetry, the Tutte polynomial 
is also a polynomial of the variable z. Chen and al. [2] propose the following 
intriguing conjecture: 

Conjecture 1 For finite planar self-dual lattices and for square lattice with free 
or periodic boundary conditions in the thermodynamic limit, the Potts partition 
zeros in the Re{v) > half plane are located on the unit circle \v\ = l. 

Our first idea was to prove this conjecture on the finite case on the large 
class of self dual graphs as possible. But later we realize that this conjecture is 
not true in the finite case. Let us take the following simple example: choose the 
self dual graph C with four vertices and six edges defined as a cycle of length 
four and an edge of multiplicity three, see Figure [TJ The Tutte polynomial of 
this graph on the previous hyperbola with q = 16 can be written as a function 
of z as 

T(C, z) = z 3 + 6z 2 + 12z + 218. 




Figure 1: Graph C : a cycle with an edge of multiplicity three 

There are one not positive real root (—2 — (210) 1 / 3 ) and two conjugated 
roots with not negative real part ((210) 1/3 /2 - 2 ± z(210) 1/3 V3/2); then, from 
the relation between z and i>, the complex roots in z with Re(z) > correspond 
with the roots in v in the region Re(v) > not on the unit circle. At least in 
the finite case, this conjecture is false. 
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In this work, at the thermodynamic limit, the analyticity of the pressure 
can be obtain by studying and avoiding the location of the accumulation sets 
of zeros of the partition function (or Tutte polynomial) for a family of self-dual 
graphs. 

A central role in our work is played by a theorem on analytic functions due 
to Beraha, Kahane and Weiss. Let D be a domain (connected open set) in the 
complex plane, and let ct\, . . . , olm, Ai, . . . , Am be analytic functions on D, none 
of which is identically zero. For each integer n > 0, let define 

M 

fn(z) = J2^k(z)[X k (z)] n . 
fc=l 

We are interested in the zero sets Z(f n ) = {zefl: fn(z) = 0} and in particular 
in their limit sets as n — > oo: 

- liminf Z(f n ) = {z e D : every neighborhood U 9 z has a nonempty intersec- 
tion with all but finitely many of the sets Z(f n )}. 

- limsup Z(f n ) = {z e D : every neighborhood U 9 z has a nonempty intersec- 
tion with infinitely many of the sets Z(f n )}. 

Let k be a dominant subscript z if |Afc(z)| > |Aj(z)| for all I e {1...M}. 
Then the limiting zero sets can be completely characterized as follows: 

Theorem 2.1 Let D be a domain in C, and let Ai, . . . Am, ai, ■ • ■ , aji be ana- 
lytic functions on D , none of which is identically zero. Let us further assume 
a no-degenerate-dominance condition: there do not exist subscripts k ^ k such 
that A fe = uj\ k > for some constant u) with \u\ = 1 and such that {z G D : 
kis dominant atz} (= {z e D : k! is dominant atz}) has nonempty interior. 
For each integer n > 0, define f n by 

M 

fn(z) = J2<Xk(z)\ k (z) n 
fe=l 

Then liminf Z{f n ) = limsup Z(f n ), and a point z lies in this set if and only 
if either 

(a) There is a unique dominant subscript k at z, and a k {z) = 0; or 

(b) There are two or more dominant subscripts at z. 

Note that case (a) consists of isolated points in D, while case (b) consists of 
curves (plus possibly isolated points where all the Afe vanish simultaneously). 

We also recall the classical famous Vitali's convergence theorem useful to 
obtain analyticity properties: 

Theorem 2.2 Let p n (z) be a sequence of functions, each regular in a region 
D. Assume that it exists a constant B as \p n {z)\ < B for every n and for all 
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z £ D. If p n {z) tends to a limit as n — > oo at a set of points having a limit 
point inside D, then p n (z) tends uniformly to a limit in any region bounded by 
a contour interior to D: the limit therefore being an analytic function of z. 

3 General results 

Here, we are interested in the following family: 

M 

f n ( Z ) [a+ cor +0 ak (z)[\- k (z)r 
fe=i 

where A+ fc (z) and A~ fc (z) are the solutions of the following equation 

X 2 - (z + 2 + a k )X + z + q + 1 = 
with a k G [0, q\. They can be expressed as follows 

A i ( z ) = \ ( z + 2 + a k ± v / (2 + afe) 2 -4(g-a fc )) ■ 

By convention, we denote A 0fc (z), the solution between A+_(z) and X~ (z) with 
the greatest magnitude. The functions {a ak ,f3 ak , k = 1..M} are chosen such 
that the function f n stays a polynomial function in the variable z. We denote 
by a u — sup fe=1 M a k and ai — vaik=i...M £Jfe- We also use some subsets of the 
complex plane: 

^]-oo,-o„-2[ = {z = c + id; (c, d) G TR 2 , c < -a u - 2}, 

D]-a,,+oo[ = {z = c + id; (c,d) G H 2 , c > -a;}. 
Then, the analyticity of the pressure in some regions can be derived easily: 

Theorem 3.1 There exists only one dominant eigenvalue at z: 

- \fz G £>]_ a ,, + oG[ \{c£ [-a/, sup(-a/, -a u + 2^q - a u )], d = 0}, A„ u (z) is the 
dominant eigenvalue. 

- Vz G -D]-oo,-a„-2[\{c G [inf(-a u -2, -a ; -2^9 - ai), -a„-2], d = 0}, A a , (z) 
is i/ie dominant eigenvalue. 

Proof: It is a direct consequence of lemma 16.11 and 16.21 

Now taking p n (z) = — JjlLJl we can deduce the following corollary. 

n 

Corollary 3.1 We have 

• //Vz G D]_ au+oo[ \{c G [-<z;,sup(-a;, -a u +2^/q - a u )], d = 0}, a au (z) ^ 
then Pn{z) — > ln[A Qu (z)] as n — > oo. 
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• IfVz e £ ) ]-oo,-a„-2[\{c e [inf(-a„-2, -ai~2y/q - a/), -a u -2], d = 0}, 
Pa t (z) ^ then p n (z) — ► ln[A ai (z))] as n — > oo. 

Let K be a non negative real and -B((0, 0), K) the disk center at (0, 0) with ra- 
dius K . Both limits are analytic functions of z respectively on subsets -D]- a ,,+oo[\ 
{ce [-a;,sup(-a;, -a u + - a u )], d = 0} f| -B((0, 0), K) and D]_ 00 ,_ au _ 2 [ \ 
{ce [inf(-o u -2,-oi-2Vg^50,-a tt -2], d = 0} Q #((0, 0), if). 

Proof: From the previous theorem, Vz € Di_ aii+00 r\ {c G [— a;, sup(— a;, — a u + 
2-^/g — a u )], d = 0}, the dominant eigenvalue is A 0u (z). Because of the given 
form of \ au (z) we obtain that 



with Ci(K) is a constant depending only of K,a u and We have also that 
Mz £ £ , ]- 00 ,-a„-2[ \ {c € [inf(-a„ - 2, -a; - 2-^/g - ai),-a u — 2], d = 0}, the 
dominant eigenvalue is A Qj . We find in this region that 



with C2{K) is a constant depending only of K,ai and q. Now the result is a 
direct consequence of Vitali's convergence theorem. 

4 Examples 

After working in two specified regions of the complex z plane, it may be be in- 
teresting to know what happens elsewhere. It is the goal of this section. Several 
simple examples of self-dual graphs are given: the classical one is the wheel. 
The other graphs - the strip of triangles with a double edge, the cycle with 
one multiple edge - are built from the wheel by moving some of its edges. By 
making this transformation, we only have to keep the symmetry property of the 
Tutte polynomial. We provide their Tutte polynomial using the contraction and 
deletion rules useful to obtain recurrence formula. We need these recurrence for- 
mula to identify without ambiguity the functions A^ fc introduced before. These 
graphs belong to the framework we discussed before for a particular choice of 
the parameters at like 0, 1 or q. We present the location of accumulation sets 
of zeros and more precisely the curves describing the degeneration of the dom- 
inant eigenvalue using respectively the variable z and v in the whole complex 
plane. For all graphs G, the notation T(G) for the Tutte polynomial T(G, x, y) 
evaluated at the point (x, y) shall be used. 

4.1 Triangles with a double edge 

These graphs are called G n (see Figure [2]) and graphs having n triangles Tr n . 
We have these relations: 



T(Tr n ) = xT(Trn-i) + T(G n -i) \ yT{Tr n ) = xyT(Tr n -i) + yT{G n . x ) 



1 < |A„ U (*)| <C X {K) 



a u - ai 
2 



< |A ,(z)| <C 2 (K) 
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Figure 2: G n a strip of triangles with a double edge 



This implies that 

T(G n ) = (1 + y)T(Tr n ) - xyT{Tr n -x). (1) 

We prove easily the following recurrence formula for Tutte polynomial of graphs 
Tr 

T(Tr n ) = {l + x + y)T{Tr n -x) - xyT{Tr n - 2 ). (2) 
Using ((T|) and it comes 

T(G„) = (1 + x + y)T(G„_i) - xyT(G n _ 2 ). (3) 

We deduce that: 

/ T(G n ) \_ { l + x + y -xy\fT(G n - 1 
Un ~{ T(G n -i) J { 1 T(G„_ 2 

The eigenvalues of this matrix are 



x + y + 1 + y/TL + x + y) 2 — 4xy x + y + 1 — Wji + x + y) 2 — Axy 
Mi= g W = 2 ■ 

The matrix defining eigenvectors is denoted by P and its inverse P^ 1 . They are 
given by : 

p=( 7 7 V ^ x = ^— f \ 

Now, it comes 
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where 

U = ( Ml + Ml + MlM2 - Ml - M2 

Then, the Tutte polynomial associated with G n is 

Mi - M2 M2 - Mi 

On the hyperbola (x — l)(y— 1) = g, these eigenvalues are of the form introduced 
before with a = 1 and the Tutte polynomial can be written as 

T(G n ) = ai (z)[X+(z)] n+1 +p 1 (z){\-(z)} n + 1 
With a ^ z) = xtW-KW ^ 01{Z) = ATW-xj-W - 

It is outstanding that such a simple form of Tutte polynomial leads to a large 
variety of cases according to the values of q when studying the degeneration of 
the dominant eigenvalue. More precisely, we have both following propositions 
respectively using complex variable z and v. 

Proposition 4.1 For the family of graphs (G n ) n >o, the location of the degener- 
ation of the dominant eigenvalue is described in the complex plane using variable 
z = c + id as follows: 
- Vq > 2 

\X+(z)\ = \\^(z)\^d = Q,ce [-1 - 2^/q~l, -1 + 2y/q~l] 
-V?G [1,2] 



lA+coi = \K{z)\ 



d = 0, c € [-1 - 2^q~l, -1 + 2^q~~T\ 
orze C((-q- l,0),2-g) 



where C((—q — 1, 0), 2 — q) denoted the circle of center (— q — 1, 0) and of radius 
2-q. 



Proof: we just have to apply lemma I6T2I in the particular case a = 1. 
By using lemma 16.41 and lemma 16.51 it comes 

Proposition 4.2 For the family of graphs (G„)„>o, the location of the degen- 
eration of the dominant eigenvalue is described in the complex plane using the 
variable v = re as follows: 
We have |A^(z)| = |A^(z)| when: 



ifq€ [1,25/16] 



9 e [arccos( 2 ^3 ),arccos( 2 j£ )], r 

or v eF(C((-q- 1,0), 2 - q)). 
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9 6 [arccos( 1 ')> 7T \j f = 1 

or 8 = ir, r G [1/ri, ri] 
or«eF(C(H-l,0),2- g )), 

S [arccos( '^VP 1 ), tt] , r = 1 
or 9 = tt, r G [1/n; ^i] 
wi/i ri and 1/ri i/ie roots o/ i/ie polynomial ^/qr 2 — (1 + 2\/g — l)r + ^/g = 0. 

The Figure [3] shows for this family of graphs, the location of the degeneration 
of the dominant eigenvalue in the v complex plane for different values of q. We 
can notice the particular value q = 1 + Vi/2. 





4.2 Wheel graphs 

We consider the wheel, a planar graph with n + 1 vertices with a cycle of length 
n and one internal vertex linked to all the vertices of the cycle (see Figure 2]). 
Then, we obtain a wheel graph having n triangles, denoted B n . The Tutte 
polynomial for these graphs can be deduced from the relations between the 
Tutte polynomials of G n and Tr„. We have 

T{B n ) = r(Tr n _x) + T(B n _ x ) + 2/(T(G„_ 2 ) - xT(G n . 3 )) 
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Figure 4: The wheel 



and 



T{B n ) - T(B n -!) = T(G„_i) - yT(G„_ 2 ) + y(T(G„_ 2 ) - zT(G n _ 3 )) 

= T(G?„_i) - xyT(G n . a ). 
The recurrence formula (|3|) provides that 

T{B n ) = (2 + x + y)T(B n ^) - (1 + x + y + xy)T(B n _ 2 ) + xyT(B n _ 3 ) 

thus 



Mn > 3, 
The eigenvalues are 



T(B n ) 
T(B„_i) 
T(B n _ a ) 



2 + x + y _(l + a; + i/ + icy) xy \ / T(B n _i) 
1 T(B„_ 2 ) 

1 / I T(B n _ a ) 



+V+1 + J (1+x+y) 2 -ixy x+y+l-J (1+x+y) 2 -4xy n 
PJ. — 2 ' ^ 2 — 2 ' ^ 3 = ^ 

The matrix defining eigenvectors P and its inverse P _1 are given by: 



P = 



Mi 



/'2 



and 



(Mi 


- l)(A*a - 


1) 


(Mi 


- i)(a*i - 

Mi 


M2) (M2 


-1)0*2 
M2 


- Mi) 


(Mi 


- l)(A*a - 


1) 


(Mi 


-l)(Mi- 


M2) (M2 


" 1)^2 


- Mi) 


V (Mi 


- 1)(M2 - 


1) 


(Mi 


-l)(Mi- 


M2) (M2 


-1)(M2 


-Mi) / 




P- 1 


-( 


' 1 


-Mi - M2 
-M2 - 1 
-Mi - 1 


MlM2 \ 
MlM2 
MlM2 / 
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Then, it comes 

/ T(B n ) \ ( 1 \ / T(B 3 ) \ 
Vn > 4, T(B„_i) = P \P ' 1 T(B 2 ) 

V t(b„_ 2 ) y v o o M r 3 / V ?W / 

where 

T(B 3 ) = x 3 + y 3 + 3x 2 + 3y 2 + 4xy + 2x + 2y 
T(B 2 ) = x 2 + y 2 + xy + x + y 
T(Bi) = xy 

We can conclude that 

T(B n ) = (xy-x-y-l) + tf+^. 

On the hyperbola (x—l)(y—l) = q, these eigenvalues are of the form introduced 
before with a = 1 for {^1,^2} an d a = q for the eigenvalue 1. The Tutte 
polynomial can be written as : 

T(B n ) = {q- 2)[\-{z)r + [A+(*)] n + {K(zT 
where A~(z) = 1. 

Proposition 4.3 For the family of graphs (B n ) n >o, the location of the degen- 
eration of the dominant eigenvalue is described in the complex plane using the 
variable z as follows: 

- Vq e [1,5], q + 2, when |A+(z)| = \X{(z)\ > 1 or |Ai(z)| = 1 

{c <= [-g, -1 + 2Vg __ T] and d = 
or 
ce[-(« + 5)/2,-g] and d 2 = -(c + <?) 2 |±§±f 

- yq > 5, w/ien |A+(z)| = \\i(z)\ > 1 

=>c£[-l- 2yjq - 1, -1 + 2yfq - 1] andrf = 0. 
-Forq = 2, when \\\{z)\ = \\{{z)\ 

4ce [-3,1], d = 0. 
Proof: Let us notice that if we denote 

Ai(*) = 1 + \f~q--\pe ia 

with p > 1 then 

|Ai(z)| = 1 + 2y/q- lpcosa + (g - l)p 2 . 
If |Ai(z)| = 1, it means 

cos(a) = — — lp/2- 
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As p > 1, the last equation has no solution for q > 5 which explains the three 
cases. 

For the second and third sentences, we apply directly the Lemma 6.2. For the 
first one, we have two different cases of degeneration. By using the Lemma 6.2, 
\Xf (z)\ = \Xi(z)\ gives the set c = [-1 - 2y/q~T, -1 + 2y/f^l\ and d = 0. 
But, on the interval [—1 — 2y/q — 1, —q], the eigenvalue 1 is dominant and does 
not match with a degeneration case. 

In the same way, the Lemma 6.3 leads us to consider the case |Ai(z)| = |A 9 (2;)|. 
On the set c e [-(q + 5)/2, -q] and d 2 = -(c + q) 2 llXl+l ' A «( z ) = 1 and can bc 
dominant with the eigenvalue Ai(z). However, on the set c G [— q, —(q + l)/2] 
and d 2 = — (c + ?) 2 ffq^pf , A g (z) = z + q + 1 and can be dominant with the 
eigenvalue \i(z) but not appears in the Tutte polynomial. In this case, Xi(z) is 
the one and only dominant eigenvalue. This is not a case of degeneration. 




Figure 5: The limiting zero sets for different values of q for the wheel in the v 
complex plane. 

Proposition 4.4 For the family of graphs (B n ) n >a, the location of the degen- 
eration of the dominant eigenvalue is described in the complex plane using the 
variable v = re %9 as follows: 

if a Gil 4] a + 2 { ^ ^ [arccos(^^), arccos(-^)], r = 1 
1 ' J ' T ' \ orv e F{T). 
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e [arccos(^4^p±),7r], r = 1 

or0 = n,re[ ^-f = * , ^ + f=* ] 
orveF(V). 

6 e [arccos(^^^±), n], r = 1 

or 9 = 7r, r G [1/ri; ri] 



€ [arccos(^=),7r], r = 1 



-For = 2<^ . 2 Yl ATI 

| or 6» = 7T, r e [^; V2] 

where ri, l/r*i are the roots of the polynomial ^Jqr 2 — (1+2^/ q — l)r+^/q = and 
= {(c,d) elR 2 ,ce [_(g + 5)/2,-g], d 2 = -(c + a) 2 (2c + q + 5)/(2c + a+ 1)}. 



Proof: We use the proof of the lemma [6741 to determinate the transformation of 
the segment [—q, — 1 + 2-^g — 1J and this one of [—1 — 2y/q — 1, — 1 + 2y/q — 1J. 
Next, for the set P, we have only modified the value of d to apply the lemma RT5l 



In the Figure [5j we show the limiting zero sets for different values of the pa- 
rameter q. The dotted line represents the degeneration case |A]*~(z)| = |A^(z)| > 
1. The other curve concerns the case where the eigenvalue 1 is dominant. When 
q is greater than 5, the eigenvalue 1 is never dominant. 



4.3 Cycle with one multiple edge 

We consider graphs having 2n edges. We assume that, for such a graph, we have 
one cycle G n+ i of length n + 1 and one edge with order of multiplicity k = n. 
We denote these graphs by G n ^. 

Let us introduce the following graphs A nt k, B n ,k an d C n ,k'- A n ^ is the graph 
having one cycle of length n and one vertex having k loops; B n ,k is the graph 
having one isthmus on length n and ended by k loops; and C n ,k is the graph 
having one isthmus of length n + 1 and the last edge with order of multiplicity 
k. These graphs are presented in Figure [Sj 
We have 

T(G n ,k) — T(C n _i,fc) + T(G n -i,k) 

= T(G n _i i fc_i) + T(_B„_i ! fe_i) + T(G n -i,k) 

= xT{C n . 2 ,k-i) + T(B n _i ifc _i) + T(G„_i, fc ) 

= a;[T(G n _i, fc _i) - T(G n - 2 ,k-i)\ + T^-i.fc-i) + T(G„_ 1>fc ) 

= x[T(G n -i : k-i) - T(G n - 2 ,k-2 - T(A n ^ 2 ,k-2)} + 

+ T(B n _i^ k _i) + T(G n -i t k) 

then 

T(G nt k) = x[r(G„_i, fe _i) - T(G n _ 2: fc-2] - xT(A n _ 2 ,k-2) 

+xyT{B n _ 2 ,k-2) + T(G n _i )fc _i) + T(A n _ hk _ 1 ) 

We deduce that 

T(G nt k) — T(G n ^i t k--L) = x[T(G n -i t k-i) — T(G„_2,fc-2] 

+ (y - x)T(A n _ 2 , k _ 2 ) + y(x + l)T(B n _ 2jfc _ 2 ). 
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Fi gure 6: From left to right and from top to bottom, the graphs G n) k, A n k, B n ,u and O n fc 

Moreover, we have both following relations 

T(A n -i i k-i) = yT{A n - 2 , k - 2 ) + ?/T( J B n _ 2 , fe _ 2 ) (4) 

T(B n -i,k-i) = xyT(B n -2.k-2)- (5) 
Taking U n ,k = T(G, hk ) - T(G n _ hk _i), it comes that 

U n ,k = xU n -i ik -l + (y~ x)T(A n -2M-2) + V{X + l)T(S n _ 2 ,fe-2) (6) 

U n -i,k-i = xU n - 2 ,k-2 + (y- x)T(A„_3 ifc „ 3 ) + y(x + l)T(S„_3 ife _ 3 ) (7) 

Using the relation and computing ©-J/©, we obtain : 

U n ,k ~(x + 2/)[/ n _ lifc _! + xyU n _2,k-2 = (xy - l)T(B„_ 2 ,fc-2)- 

Considering this last relation with order n — 1, k — 1, the equation ((5|) and after 
subtracting and multiplying by xy, it leads to 

U n ,k - (x + y + xy)U n -i ik -i + xy(l + x + y)U n - 2 ,k-2 - (xy) 2 U n - 3ik - 3 = 0. 

Now, writing with the help of previous Tutte polynomials, we find the following 
relation: 

T(G„,n) = TT(G n - lin .. 1 )-(rf-l + xy(l + x + y))T{G n -2,n-2) 

+JxyT(G n - 3 . n ^ 3 ) - {xy) 2 T(G n ^ i . n ^ 4 ) 
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where "f — xy + x + y+ l. Hence 

/ 7 
1 


V o 



T(G n -2,n-2) 



(7 — 1 + xy(l + x + y)) jxy 



1 
1 



(xy) 2 \ ( T{G n -\ lTl -\) \ 



J 



T{G r , 
T{G r , 
\ T(G r , 



1,71 

2, n-2 

3, n — 3 

4, n-4) ) 



( T(G n 


n) 


> \ 




( 1 








^ 


T(G n 


-1,71- 


-l) 


= P 





X 








T(G n 


-2,71- 


-2) 








y 





V r(G„ 


-3,n- 


-3) J 




V 








xy ) 



71-3 



p- 



/ r(G 3 , 3 ) \ 

T(G 2 . 2 ) 

V r(G ,o) / 



where 



T(G 3 , 3 ) = x 3 + y 3 + xy 2 + a; 2 y + xV + a; 2 + y 2 + xy + x + y, 



T(G 2 , 2 )=x 2 +y 2 +xy + x + y, T(G M )=x + y, T(G ,o) = 1. 
Then, the Tutte polynomial is given by 

xy-x-y (xy) n 



T(G n 



(x-l)(y-l)' 



(x-l)(y-l)' 



On the hyperbola (a; — l)(y — 1) = q, these eigenvalues are of the form 
introduced before with a = for {x, y} and with a = q for {1, xy}. The Tutte 
polynomial can be written as : 



T(G„,„) = i[A+(z)r + i-^[A-(z)]" 



'Aj(z)]' 



9-1 



Proposition 4.5 For i/ie family of graphs (G„ jn ) n >o, fie location of the de- 
generation of the dominant eigenvalue is described in the complex plane using 
respectively variable z and v as follows: 
\fq > 1, when \X (z)\ = |A,(z)| 



(c, d) G iR 2 , c G [-2 - g/2, -g/2] , d 2 

«g A_^ /2 uC((-i/Vg,o),i/Vg) 



-(2c + g + 4) 



(c+g) 



2 1 



2c + g 



where A_ y ^/ 2 *s fie line Re(v) — — y/q/2 and C((— 0), denotes the 

circle of center (— l/«/g, 0) and 0/ radius 1/ ^/g. 



Proof 1 TTie /irs£ pari follows again from lemma flOl with a — 0. ./Vow, using 
variable v, it is easy to see that \y\ = \l+^/qv\ = 1 Zearfs tovE C((— 1/^/g, 0), 
and \x\ = |1 + y^g/ul = 1 to v <E A_ v ^/ 2 . 

We notice that, for g > 4 or — g < — 2^/g, then the unit circle |u| = 1 is not 
at all in the set of degeneration of the dominant eigenvalue. This is completely 
the opposite of what we are waiting for in the case of self dual strip of the finite 
square lattice. 
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5 Concluding remarks 



1. Given a Tutte polynomial written with M pairs of eigenvalues like (A+ fc {z), X~ k (z)), 
k = X...M. The cases of degeneracy look like to |A„ (z)| = |A~ (z)| or 
|A+ {z)\ = |A~,(z)|, for some k,k' in {1,...,M}, where these eigenvalues are 
dominant. As the products |A+ fc (z)X~ k (z)\ are constant, the first case implies 
that all eigenvalues have the same magnitude. By using the lemma 16.21 the 
degeneracy region in z is given by C\ [— — 2^/q — Ok~; —a-k + 2y/Q — &k\- 

k=l...M 

For the second case, by using the proof of the lemma UTTl |A+(z)| is a not in- 
creasing function of a, a £ [a; ;a u [, on the set Lh_ oc ._ a i and a not decreasing 
function of a, a € [a/ ;a u [, on the set .D]_ ai . +00 [. Thus, we have to study only 
one case of degeneracy of these sets : |A ai (z)| = |A au (z)|. 

2. To conclude, in the finite case, the conjecture of [5] is true for the strip of tri- 
angle with a double edge and the wheel but it is false for the cycle with an edge 
of high order of multiplicity. But, as a direct consequence of theorem 13.11 we 
conclude that, for all graphs studied of this paper, the accumulation set of zeros 
in positive half plane Re(v) > is only a sector of the unit circle |i>| = 1. Notice 
that from corollary 13. 1[ we even derive an analyticity result on the positive real 
axis as it is foreseeable for some kind of graphs. But in the general case without 
knowing other explicit form of eigenvalue, a complete description of the set of 
degeneration of the dominant eigenvalue in the positive half plane remains open. 

3. In the other hand, we found many other self dual graphs belonging to this 
framework: these graphs can be seen as a transition between the wheel and 
the cycle with a multiple edge. Moreover, we point out that other choice of 
parameter a might be interesting. For example consider a family of self dual 
strip graphs defined as follows. Let the L x x L y lattice strip have periodic 
boundary longitudinal or horizontal boundary condition and connect all of the 
vertices on the upper side of the strip to a single external vertex, while all of the 
vertices on the lower side of the strip have a free boundary condition. The case 
L y = 1 is the wheel graph discussed previously. For this family of graphs, the 
Tutte polynomial is already known, in particular for the last but one transfer 
matrices, the eigenvalues obtained for L y = 2 and L y = 3 are respectively roots 
of the following equations 

A 4 - (20 + y) + 3)A 3 + (x 2 + y 2 + 3(x + y) + 4xy + 1)A 2 

-xy(2(x + y) + 3)\+(xy) 2 = 

and 

A 6 - (3(:c + y) + 5)A 5 + (3(x 2 + y 2 ) + 10(a: + y)+ 9xy + 6)A 4 
-(x 3 + y 3 + 5(x 2 + y 2 ) + 9xy(x + y) + 20xy + 6(x + y))X 3 
+xy(3(x 2 + y 2 ) + 10 (a; + y) + 9xy + 6)A 2 - x 2 y 2 (3{x + y) + 5)A + {xyf = 0. 
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For this kind of graphs, the complex-temperature phase diagrams have been 
already computed for example in pQ. For these equations, we find eigenvalues 
of the form we have studied in this paper. Denoting by b r = 2 + 2cos(27r/r), the 
Tutte Beraha numbers, we just have to choose values of a = 65, 3 — 65 for L y = 2 
and a = 67, 2 — 3 — 67 + y/bj for L y = 3. Unfortunately, it seems that the 
eigenvalues of others transfer matrices was not of the desired form. It would be 
interesting to find what choice of parameters a,a a ,f3 a allows to include other 
known family of self dual graphs. 

6 Proofs 

In this section, we provide technical lemma on the degeneration of the dominant 
eigenvalue. 

6.1 Study of the magnitude of the eigenvalues 

For a given a, the eigenvalues A J" (z) and A~(z) are the solutions of the following 
equation : 

X 2 - (z + a + 2)X + z + q + 1 = 

with z = c + i d a complex not real (d ^ 0). We will discuss later on the case 
d = 0. We find 

^(*) = \ ([a + c + 2 + id±^J{z + a + 2) 2 - 4(z + g + 1)) . 
By introducing 

A = Re((z + a + 2) 2 - 4(z + q + 1)) = (a + c) 2 - d 2 - 4 {q - a) 
B = |0 + a + 2) 2 - 40 + q + 1)| = ^A 2 +4d 2 (a + c) 2 , 

the eigenvalues can be written under this form: 

A±0) = !(Wc + 2± v /i±i) + i ( d± y^( a + c) ) . 

Moreover, by using the relation B 2 — A 2 — Ad 2 (a + c) 2 , it comes 

A +o) = i u+ c +2+(-ir 
< \ 

A"0) = | la + c+2-(-iy 

where 

_ ( if a + c> 
1 1 otherwise 



f if d > 
and p = < , 

[ 1 otherwise 
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We consider the following sets 



= {z = c + i d, (c, d) e H 2 , u < c < v, d > 0} 
D ]u.v[ = {z = c + id, (c,d) e IR 2 , u < c < v, d < 0} 



We can state a lemma on the variation of the magnitude of these eigenvalues: 

Lemma 6.1 given a Tutte polynomial written with M pairs of eigenvalues like 
(\+ k (z),\- k (z)),k = l...M. ' 

Va e [a;; a u ], |A.(z)| = |A+(z)| is a not decreasing function of the variable a for 
all zeD'_ ai +oc[ for • = +,-. 

Va e [a;; a u ], \X.(z)\ = |A+(z)| is a not increasing function of a for all z e 

£f-oo,-a„-2[/°r» = +>-- 

Proof: The magnitude of A+(z) is given by: 



4|Aj(z)| 5 



a + c+2 + (-iy 



A + B \ 
2 ) 



d+(-ir 



B-A 



Taking U = \J M f 6 - and V = \J^-, it comes that 

(2|A+(z)| 2 )' = (a + c + 2 + (-l) n V)(l + {-l) n V) + (d + {-lfU){-lfU' 
then 



(2|A^(z)| 2 )' = (a+c+2+(-l)"y) (l + (-i r ^±^y { -i nd +(-irU) 
Moreover B' = ±{AA' + 4d 2 (a + c)) = -f . We find that 



B'-A' 



(8|A a + (z)| 2 )' 



VB 



[a + c + 2 + (-l) n V] [AVB + (-l) n (A'B + W)] 
+ { - T ^\d+(-iy U ](W-A'B). 



Next, we express this differentiation in function of U and V by computing: 

' \a + c\ = UV/\d\ thus a + c= (-l)"+*W/d 
B = U 2 + V 2 
A = V 2 -U 2 



A'B + W = AV 
W - A'B = AU 



(-l) n +"UV 2 
d 

( (-l) n+p U 2 V 

{ d 



+ 2V + d(-l) n+p U S j 
2U + dV(-l) n+p 



We conclude that 
(2|A^(z)| 2 ) 



2y_ (-1)"^ 



B 



VI 



U(-1)p 
d 



[u + i-ifdy 
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Hence it is enough to remark that not increasing or not decreasing depends only 
on the parity of n. Now, we identify the eigenvalue between A+(z) and X~(z) 
with greatest modulus which will be denoted by X a (z). We have 



\>i(z)\=\X-(z)\=> (_l)n (a + c+2 )^±£ = -d(-l)P ' B A 



Thus, by looking at the sign of the members of the previous equalities: 

\Xt(z)\ = \X-(z)\^\ Z z l^+ a - 2 '- al 
It comes by the monotonicity of the eigenvalues 



X+(z)onD^_ oo _ a _ 2[ 
Xi(z)onD^ a+oc[ 



X a (z) = 

Given a Tutte polynomial like 

M 

f n (z) = £[a 0fc (*)[A+ (*)]" +0a k (z)[X- k (zT 
fc=l 

and ai = min fe= i.. M a k , a u = max fe= i.. M a k - 

We can assert that, for ai < a < a u , |A.(z)| is a not decreasing function of 
the variable a for all z G D'_ a[ +oo j = D'_ ak +oc j. But also that, for 

k=l...M 

ai < a < a u , |A.(z)| is a not increasing function of a for all z G D' ^ _ 0u _ 2 [ = 
P D"_ 0O| _ Ofc _ 2 [ for • = +, — . This conclude the proof. 

k=l...M 

Remarks: 1) The case d = implies that, if (a + c) 2 — A(q — a) > 0, the 
eigenvalues are reals and differents. We can not have |A+(z)| = |A~(z)|. How- 
ever, if (a + c) 2 — 4(q — a) < 0, the eigenvalues are conjugated and we have 
|A a (z)| = lA^z)! = |A+(z)|. This gives a case of degeneration and it will be 
discussed in next lemma. 



2) We notice that, on the set D'_ a 2 _ o j, the eigenvalue X a (z) can be 
identified to be A+(z) or X~{z). So, we can have a degeneration case like 
| A~ (2) I = |A+(z)| or, when we have M couples of eigenvalues, a degeneration 
case like |A+ (z)\ = \X~. (z)\ for i ^ j. These cases are studied in the next lemma. 

6.2 Study of some degeneration cases 

At first, we consider a function using one and only one pair of eigenvalues. 
In this case, it should be interested to study the case of the equality of the 
magnitudes: 
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Lemma 6.2 Va e [0, q - 1] 

|A+(^)| - \\~(z)\ ^d = 0,ce[-a- 2^9=5, -« + 2^~^ 
Va G [a - l,q] 



\X+(z)\ = \X-(z)\ 



d = 0, c € [—a — 2^/g — a, —a + 2-^/g — a] 
orze C((-g - 1, 0), 1 - q + a) 



where C({—q — 1, 0), 1 — q + a) denoted the circle of center (— g — 1, 0) and of 
radius I — q + a. 

Proof: The equality in magnitude is true only for z € -D*_ a _ 2 _ a [ for • = +,—. 
We have to solve this equation 



(-!)"(« + c + 2 ) V ^ = V 

It conies 

A[(a + c + 2) 2 + d 2 } = B[d 2 - (a + c+2) 2 ]. (8) 
As these terms must have the same sign, we hnd 

A 2 [{a + c + 2) 2 + d 2 } 2 = [A 2 + 4d 2 (a + c) 2 ] [d 2 -{a + c + 2) 2 ] 2 

thus 

d 2 [A(a+c+2)-{a+c)(d 2 -{a+c+2) 2 )][A(a+c+2) + {a+c)(d 2 - (a+c+2) 2 )] = 
It can be written as 

d 2 [A(a + c + 2) - (a + c)(d 2 - (a + c + 2) 2 )] [(c + a + l) 2 + d 2 - {q - 1 - a) 2 ] = 

The case d = implies that the eigenvalues must be conjugated and gives the 
segment c G [—a — 2^Jq — a, —a + 2y/q — a\. 

The third term of the product provides the equation of the circle C((—q — 
1, 0), |1 — q + a\). But, as —a — 2 < c < —a, this makes sense only if 1 — q + a > 
thus q — 1 < a < q. 
The second term implies that 

d 2 = A(a + c+2)/(a + c) + (a + c + 2) 2 . 

If we put it into the equation ([5]) , it follows that 

^[(a + c+2) 2 + d 2 ] = AB(a + c + 2)/(a + c). 

As —a— 2 < c < —a in the set D?_ a _ 2 _ Q r, for this value of d, both previous terms 
have opposite sign, then it is not a solution. That ended the proof of the lemma. 
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Remark : An other way to prove this lemma is to notice that 



{ 



|A+(*)A-(*)| = |* + g + l| 
(X+(z)-l)(X-(z)-l) = q-a. 



Thus, it is natural to consider the following form of the eigenvalues: 



A+(*) 



1 + V? - ape la \ a (z) = 1 + V? - a/pe 



We obtain 



|A+(*)I = |A-(z)| ocos(a) = -^O(p+i) rp = 1. 

z p 



The case p = 1 describe the situation where the eigenvalues are conjugated. 
This means that 



In this remark, we have just proved the first implication of the previous lemma. 
It turns out that these expressions of the eigenvalues are useful to study the set 
of zeros for our family of finite graphs and gives directly an idea of the location 
of zeros. 

Now, consider a function using at least one pair of eigenvalues like (A+, A~) 
and the pair (A+ , A~ ) . The degeneration can happen when the magnitudes of 
one eigenvalue of each pair are equal: 

Lemma 6.3 



Proof: First, we remark that \{z) = 1 or \{z) = z + q+1. Thus, the relation 
|A (z)| = |A g (z)| becomes 

|A a (z)| = 1 > \z + q+ l|or|A a (z)| = \z + q+l\ > 1. 

In all cases, we have to study |A+(z)| = 1 or |A~(z)| = 1. It comes 



d = 0, c e [—a — 2y/q — a, —a + 2y 'q — a}. 



The other case is possible only when a > q — 1. 
We deduce that if p ^ 1, for these value of cosinus, 





X t( z ) + K( z ) = c + 2 + a + id^ A+(z) + l + 



q — a 



c + 2 + a + id. 



A+(z)-l 
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Taking, for example, A+(z) = e la , 
q — a 



e la + 1 



2(1 — cos a) 



(cos a — 1 — i sin a) = c + 2 + a + id. 



It gives 



c = cos a — 1 — g/2 — a/2 
g — a 



d = sin a 



1 



2 — 2 cos a 



As cosa = c + g/2 + a/2 + 1, assuming that — (q + 4 + a)/2 < c < — (g + a)/2, 
we deduce that 



d = ±\/ 1 — cos 2 a 



1 + 



q — a 



2c + q + a 



-(c + q) 



2 2c + g + 4 + a 
2c + g + a 



6.3 Correspondence between z and v, case d = 

Many authors work with the variable t> = re such as z = c + id = ^/q(v + 1 / v) . 
It leads to c = ^(r + 1/r) cos(#) and d = ^(r — 1/r) sin(#). We call F the 
transformation from zto ti: 



v = F{z) = 



± V^ 2 - 4g 



We have to calculate F(z) when z e {(c, d) e 1R 2 , d = Oandc e [— a — 
2^/g — a, — a + 2^/g — a]} and when z belongs to the circle 

C((-q - 1, 0), 1 - q + a) = {(c, d) e IR 2 , (c + q + i f + d 2 = (1 - q + a) 2 } 

for g — 1 < a < q. 



Lemma 6.4 Let define 
9(a) -- 



-9(a)+V g (a) 2 -4 _ -h(a) + ^/ h(a^ -4 

' 1 2 ' 2 2 



• Vg>l,Va< 4(^-1), 
d = 0, c e [(/(a)^, /i(a)7g] 

• VI < g < 4, Va > 4(^-1), 

d = 0, c e [.9(0)^/9, M«) V?] 



fle[arccos(^),7r] q f 6 = n 

r = 1 ( r£ [l/ri;n] 



G [arccos(-r^), arccos(^p-)] 
r = 1 
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. Vg>4,Va>4(V9-l), 

d = 0, c G [o(a),/g, /i( a )\/<7l I m / 1 / i , , r i 
un /vi' wv*j | r G [l/r 2 ; 1/ri] U [ri;r 2 ] 

Proof: Taking d = = — 1/r) cos(0) and c = y^(r + 1/r) cos(0), we must 
have 



g(a) = --2J1-- < (r + 1/r) cos(0) < /i(a) = - + 2.1--. 

V i V i 

Remark that g is not increasing on [0, q—1] (<?(0) = —2, g(q — 1) = — (^fq + 
l/y/qj) and not decreasing on [q — l,q] (g(q) = —^/q)- 

h is not increasing on [0,q] and h(0) = 2, h(4[y/q — 1]) = —2 and h(q) = —^fq. 
We hrst look for the values of 9 when r = 1. We have to solve 

g(a)/2 < cos(0) < ft(a)/2. 

- For q > 4, we have g(a)/2 < -1. Thus, if a < 4(^ - 1), G [arccos(^), tt] 
and if a > 4(y / g— 1), h(a)/2 < —1 which gives no solution. 

- For q < 4, we have g(a)/2 > -1 and g(4{y/q- 1)) = -2. Thus, if a > 4(^-1), 
6» e [arccos(^),arccos(2M)] and if a < 4(^-1), G [arccos(^), tt]. 

Now, we look for the values of r when 9 = ir. We have to solve 

-h(a) < (r + 1/r) < -5(a). 

- For g < 4, —h(a) < 2. If a > 4(^/5 — 1), —g{a) < 2 which gives no solution. If 
a < 4( v /g — 1), we obtain r e [1/ri; ri]. 

- For g > 4, if a < 4(^/q — 1), -/i(a) < 2 which gives r G [l/ri;ri]. And, if 
a > 4( v /g - 1), -/i(a) > 2 which gives r G [l/r 2 ; 1/ri] U [n; r 2 ]. 

Remark that the case = only implies r = 1. 

The following lemma gives an analytic expression of u when we apply the 
transformation f ona circle: 

Lemma 6.5 For all q — 1 < a < q, 

z = c + ide C((-q - 1,0), 1 - q + a)) ^ 

f ^(l-g + a^-Cc + g + l) 2 
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wtft A = c 2 - d 2 - 4q and B = \M 2 + 4d 2 c 2 . 
Proof: We have just to remark that, for a given z, as 

z± v / z 2 - 4q 
V = F{Z)= 2^q ■ 
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v can be written under this form 



v = 



a ± (^)-i 

2V? 




where 



B 




A 


2 


B 




A 


2 



f if c> _ f if d 

\ 1 otherwise an P \ 1 other 



> 
otherwise 
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